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A deterministic Bohmian mechanics for operators with continuous and discrete spectra is pre- 
sented. Randomness enters only through initial conditions. Operators with discrete spectra are 
incorporated into Bohmian mechanics by associating with each operator a continuous variable in 
which a finite range of the continuous variable correspond to the same discrete eigenvalue. In this 
way Bohmian mechanics can handle the creation and annihilation of particles. Examples are given 
and generalizations are discussed. 
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(N| ■ I. INTRODUCTION 

In 1951 David Bohm introduced an extension of quantum mechanics, now called Bohmian mechanics, in which 
particles have definite positions and velocities at all times, including between measurements 2 . He also showed how 
the same idea can be applied to field configurations. Bohmian mechanics agrees with all of the experimental predictions 
of the Copenhagen interpretation of quantum mechanics but has the advantages of a smooth transition to classical 
mechanics, the absence of wave-function collapse, and not having to separate the universe into quantum systems and 
classical measuring devices. Furthermore, the dynamics is deterministic and time reversible. Randomness enters only 
in the initial conditions of the particle and field configurations. On the other hand, Bohm constructed his mechanics 
f^i ■ in first quantized form which is fine for the nonrelativistic Schrodinger equation but it can not handle fermion particle 
t-H _ creation and annihilation in the second quantized representation of quantum field theory. In this paper we generalize 
Bohmian mechanics so that it can handle operators with discrete spectra, and thereby accommodate second quantized 

Ofermionic field theories while remaining deterministic. 
In Bohmian mechanics, the quantum state and Hamiltonian do not constitute a complete description of physical 
reality as they do in the Copenhagen interpretation of quantum mechanics. Bohmian mechanics singles out a set of 
Oh, dynamical variables associated with a particular set of commuting operators (for example, the positions of all of the 
i ' particles) that are needed in addition to the quantum state and Hamiltonian to constitute a complete description 
of physical reality at all times. J. S. Bell has coined the term 'beables', short for 'maybeables' for those operators 
promoted to reality status 3 . The hesitancy is built into the word to indicate, somewhat analogously to gauge freedom, 
p ; =. that the set of commuting operators is not uniquely determined by experiment. Bohm's original formulation was 
• • ■ expressed in first quantized notion His formalism required that all of the beables have continuous spectra which 
is fine when electron position operators are the beables. However, the Dirac equation presents a problem in this 
formalism since it requires an infinite number of negative energy electrons in the vacuum, so the most naive extension 
of Bohmian mechanics to the relativistic realm requires the calculation of an infinite number of trajectories even when 
there is nothing really there. If instead one thinks of the vacuum as a vacuum, and moves to a second quantized 
representation, then one is faced with the possibility of electron-positron pair creation and annihilation. This can not 
be straightforwardly treated with continuous variables alone since the number of particles is not a continuous variable. 

In 1984 J. S. Bell constructed a Bohmian-like mechanics for beables with discrete spectra, in particular he considered 
fermion configurations in relativistic field theories 3 . He gave up on particle trajectories altogether and considered in- 
stead the fermion number density operators as the beables. Unlike Bohmian mechanics, Bell's mechanics is stochastic. 
Bell was dissatisfied with this since "...the reversibility of the Schrodinger equation strongly suggests that quantum 
mechanics is not fundamentally stochastic in nature. However I suspect that the stochastic element introduced here 
goes away in some sense in the continuum limit. 3 " Since Bell there have been several other contributions to Bell-like 
dynamics for relativistic quantum field theory, some bringing back explicit particle trajectories, but like Bell regret- 
tably sacrificing determinism 4 . In this paper we introduce an alternative to Bell's mechanics that is a deterministic 
time-reversible Bohmian mechanics for operators with both discrete and continuous spectra. The determinism and 
time reversible invariance is present at the course grained level. No continuous limit is necessary, and no modification 
of the Hamiltonian need be made. Since our extension of Bohmian mechanics allows for discrete beables it is able 
to handle particle creation and annihilation in second quantized field theories, and thereby dispense with one of the 
objections to Bohmian mechanics. 

In section II of this paper we present a generalization of Bohmian mechanics for operators with continuous spectra 
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that is also amenable to operators with discrete spectra. In section III we incorporate operators with discrete spectra 
into the formalism. In section IV we present the exactly solvable case of Bohmian mechanics for one beable, and an 
intriguing visualization of Bohmian mechanics for any number of beables. In section V contains a summary of our 
results. 



II. BOHMIAN MECHANICS WITH PROJECTION OPERATORS 



In this section we present Bohmian mechanics in a generalized way, making extensive use of projection operators, 
that will allow us to incorporate operators with discrete spectra. The generalization agrees with Bohm's original 
formulation when the beables are particle position operators but also allows for any choice of commuting operators 
I = 1,2, ...L, [£,£,£,£>] — for the beables. The operators can have continuous or discrete spectra. Our generalization 
of Bohmian Mechanics is not unique. In particular, we have some freedom to choose how many and which commuting 
operators we require to describe the status of all possible measurement devices. In this section we will deal only with 
operators with continuous spectra and extend the formalism to discrete operators in the next section. 

If & has continuous spectra we can express it as 



d\tZt(\t)P t (\e) (2.1) 



where Xt parameterizes the eigenstates of the integral is taken over the entire range of Xe, and £e(Xe) is the 
eigenvalue of ^ associated with the eigenstates labeled by Xe- 

£e\\ i ,q,t>=&(\ i )\\ e ,q,£> (2.2) 

where q distinguishes states with the same eigenvalue of ^. The projection operator for the eigenstates associated 
with £e(Xe) is 

P i (X i )=^\X i ,q,e><X i ,qJ\. (2.3) 
<i 

In this expression the sum over q represents the sum or integral over states with the same eigenvalue of ^. The 
simplest &(A^) function for an operator with continuous spectra is &(Xi) = Xi in which case Xg has units of ^. For 
this case the projection operator density, Pg(Xe), takes on the particularly simple form 

P e (X e ) = S(X-ie) (2.4) 

For operators with discrete spectra, which we cover in the next section, we will find it convenient to associate a range 
of lambda variables to a single eigenvalue. 

Bohmian mechanics, describes the dynamics of a set of Xt{t)'s (that is, Ai, A2, A3..., which we denote collectively by 
A) and thereby a set of £e(Xg(t))'s (which we denote collectively by S) which represents the physical values of these 
values at time t, regardless if a measurement is made or not. The quantum probability distribution of a particular A 
configuration is conveniently written in terms of the projection operators 

L 

P(X,\t>)=<t\l[Pe(Xe)\t> (2.5) 

i=\ 

This is a probability distribution since 

J dXt J dX 2 ...j dA i /(a(Ai),6(A 2 ),...a(AL))P(A,i) =<t|/(£i,fc, ...&)!*> (2.6) 

The result is unambiguous since we require that the £ all commute with each other. The probability distribution has 
all the properties required of a classical probability distribution. The integral of the probability distribution taken 
over all of A space is 1, and the probability distribution is real and non- negative provided that all of the projectors in 
the operator product commute with each other. The projectors will commute if the ^ all commute with each other, 
which is why we made this requirement for our set of beables. 
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The quantum state is in general not an eigenstate of the S. It is propagated forwarded in time as in conventional 
quantum mechanics 

ih^f =H\t>. (2.7) 

Additionally, the A configuration is propagated forward in time with the first order equations 

d\ e (t) 



dt 



v e ({X(t)},t) (2.8) 



where the ve({X(t)},t) are chosen so that the classical probability distribution of the A configuration of an ensemble 
of identical experiments 

P c (A,t) = J dAi(O) J d\' 2 (0)...J d\' L (0)P c (\',0)J[5(X e - Aj(t,{A'(0)}) (2.9) 

agrees with the quantum probability distribution at all time provided they are in agreement at any one time. If the 
beables are chosen such that all measurements are measurements of the beables, this guarantees that the results of 
Bohmian mechanics are consistent with the results of conventional quantum mechanics. 

Using Bohmian mechanics Eq.(2.8), the time dynamics of the classical probability distribution is 

dPc({X},t) ^ dP c ({X},t)v e ({X},t) 

* " h ^ ' 

Whereas the time derivative of the quantum probability distribution on the other hand is 

^^=E < *i(n^|i[A(A/) ) fl]( n pk{x k ))\t>. (2.ii) 

i=i \j=i / \k=e+i ) 

Following David Bohm's insight, we note that if the classical and quantum probability distributions agree at any 
particular time then they agree for all time provided that the ve({X}) are chosen so that the two time derivatives Eqs. 
(2.10) and (2.11) are equal. This is what we do now. 

Our goal is to rewrite Eq.(2.11) in the form 

dP(X,t) _ f^i) 

— " dX e ' ( ' 

with Jf (X,t) real. If this can be accomplished then we can set 

«<M>- («») 

and we will have determined a consistent Bohmian dynamics. Associating the I terms in both expressions we have 

= -<t\ ^n^(Ai)j ^[Pe(Xe),H] ^ P^X^ \t > . (2.14) 
Ignoring for the moment the possibility that the right hand side of Eq.(2.14) is not real we write 

MX,t) =< t\ (J] Pt'(Xt'(t))\ MXe(t)) ( J] P A \ l >. ( 2 - 15 ) 

\i'=i ) \i"=t+i ) 



in which 



^M = -I[P (W ,H] ,2.16) 
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which is easily integrated to 



where 



MMt)) = ±r[Ge(Xe(t)),H] = ~[L t (Xe{t)),A\. (2.17) 



G e (Xt(t)) = / dX'P e (X' e ) (2.18) 

J Mi) 



and 



L e (X e (t)) = J dX'P e (X' e ). (2.19) 

are projection operators for all states greater than or less than Xg{t) respectively. The current operator, Jg, can be 
generalized to periodic beables such as the position of a bead on a ring with 

MHt)) = -^JdX' J dX"Pe(X' e )HPe(X")f(X'e, X e (t), A") (2.20) 

where f(X' e , X((t), A") = +1(— 1) if there is a non-crossing path that goes from X' e to A" through Xg(t) in the positive 
(negative) direction and f(X' t , Xe(t), A") = if there is no such path. We will not use this generalization in this paper. 

The right hand side of Eq.(2.14) is not guaranteed to be real unless [[Pi, H], Py\ = for t ^ I' which is true if 
[[&, H],&] — for I ^ £'. If this is not the case we can make it real simply by taking the real part 

MX,t) = Re ^< t\ (j[ Pt'Mt))^ Je(X e (t)) ^ [] P e ,)j [t >^J . (2.21) 

but this picks out a particular order of the operators for special treatment. A more democratic way to guarantee that 
the current is real is the symmetric average 

^(A,|t>)=<t|s|^JJA'(V(*))^^(M*))^ II A»^j|t> (2.22) 
where S{...} implies a symmetric average of all of the operators inside the braces. For example if L=3, 

Ji(Ai, A 2 , A 3 ) = I < t|2JiAA + P2J1P3 + +P3J1P2 + 2P 2 P 3 J 1 \t > (2.23) 
6 

where we have used the fact that the P's commute to combine some terms. This freedom in the choice of Jg is a 
particular example of a more general freedom. We can add any function Qe to the probability current density, Jg, 
such that ^ l d/ dXgQg — 0. This is the second way that the dynamics are not unique, (the first being the choice of 
which operators to anoint to beable status) . 
Our final expression for the A dynamics is 

dMt) =M{mht) = < t\s{(m^Mx e ,(t))) Jdxm (n^ +1 a»)} it > (2 24) 



dt 



<*IIL=iW*)|t> 



An equivalent way to write the equations of Bohmian mechanics, is to use the Heisenberg representation in which 
the quantum state does not change with time but any operator A(t) depends on time via 

A{t) = e -Ht/ l hj Le Ht/ l n_ ( 2 25 ) 

The equations of Bohmian Mechanics in the Heisenberg representation are 

<0|sj^p P ei (X e ,(t),t^J dL e (X e (t),t) JJ iV^j|0>=0 (2.26) 
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where 



in which 



dt ih 
and 

dL e (Xe(t),t) 
d\e(t) 



f , m . di e (x e (t),t) d\ t (t) di e (x t (t),t) 

dL e (X e (t), t) = + — , (2.27) 



dL t (Xt(t),t) l [Li{Xi{t) ^H] = -J e (Mt),t), (2.28) 



= P t {\t(t),t). (2.29) 



In the Schrodinger representation, since the quantum state changes with time, one is tempted to think of it as a 
dynamic variable just like the beables and wonder why the beables depend on the quantum state but not the other way 
around. In the Heisenberg representation the quantum state is not a dynamic variable so this apparent asymmetry 
does not arise. 

If an ensemble of identical experiments are performed, in which the A configurations at a particular time for each 
experiment are taken at random from the quantum distribution Eq.(2.5) then if the A configurations are propagated 
forwards and backwards in time via Eq.(2.24) or Eq.(2.26) then for all other times the probability distribution of the 
A configurations over the ensemble Eq.(2.9) will equal the quantum distribution. Provided that the set of anointed 
operators is sufficient to describe the status of all measurement devices, Bohmian mechanics will agree with all 
results of conventional quantum mechanics without resorting to wavefunction collapse or some other alternative to 
the Schrodinger dynamics to describe the measurement process as is done in the orthodox Copenhagen interpretation. 
Bohmian mechanics replaces this with the mystery of how to explain why the classical and quantum probability 
distributions should agree at any time at all 5 . 

If all of the A^'s correspond to position coordinates in nonrelativistic quantum mechanics then the conventional 
form of Bohmian mechanics is recovered from Eq.(2.24) or Eq.(2.26). For then the current operators are 

Jt(*t) = \\— ,S(X e -Xi)} + . (2.30) 
2 me 

This is the current operator in conventional Bohmian mechanics so the equivalence with the traditional formalism 
is proved. The present formalism is more flexible than the traditional formalism and can be easily generalized to 
account for beable operators with discrete spectra, which is what we consider in the next section. 



III. FITTING THE DISCRETE SQUARE PEG INTO THE CONTINUOUS ROUND HOLE 

In this section we incorporate beables derived from operators with discrete spectra into Bohmian mechanics. Doing 
so, we are immediately faced with the question of how to retain determinism which Bell sacrificed with regret. The 
problem is most apparent when the initial quantum state is an eigenstate of each ^ for then the initial S configuration 
is uniquely determined so there appears to be no room for any randomness initially. The quantum time dynamics 
will immediately make the quantum state a superposition of S states so the classical probability distribution must 
develop some spreading to agree with its quantum counterpart. Since there is no randomness in the initial conditions, 
it appears that randomness must enter through the dynamics and therefore determinism must be sacrificed. Note 
that this is only a problem if all of the £e operators have discrete spectra. If even just one of the £e operators have 
continuous spectra then that is enough to make the initial state not unique so there is the possibility of a deterministic 
dynamics producing the correct future probability distributions. Strictly speaking the same problem can arise for 
operators with continuous spectra but for continuous spectra one can wiggle one's way out by asserting that in any 
actual experiment the quantum state is never exactly in an eigenstate of all of the operators. There is always some 
spreading for whatever reason. 

A way to retain deterministic dynamics with discrete operators is to assign finite ranges of Xe to the same eigenvalue 
of & sits still when A^ is moving smoothly through a region that corresponds to the same eigenvalue and £e makes 
a sudden hop to a new eigenvalue when Xe smoothly moves from one eigenvalue region to another. The problem of 
determinism is solved by this devise since there are now many A configurations corresponding to the same physical 
state. It might appear that for operators with discrete spectra, Xe can legitimately be called a hidden variable since 
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the particular value of \g inside an eigenvalue range is unobservable and therefore " hidden" . But \? is observable 
at those times when the £e hops to a new value since there is a unique value of X( for each transition. Therefore Xe 
encodes observable information about previous transitions and the times that they occurred. We now construct the 
explicit Bohmian mechanics for the discrete case. 
If & has discrete spectra we can express it as 

ii = Y,&nPe(n) (3.1) 

n 

where is the nth eigenvalue of ^ and the projection operator is 

P e (n) = ^\n,qJ><n,q,e\ (3.2) 
<? 

, where as in the continuous case, the sum over q represents the sum or integral over states with the same eigenvalue 

of 6. 

We now seek to express Eq.(3.1) in the continuous form Eq.(2.1) and define an appropriate £,e{Xe) function and 
projection operator Pe(Xe) so that we may carry over Eq.(2.24) or Eq.(2.26) unchanged for the dynamics. There are 
many ways to do this. Technically, 6e(Xe) — Xe and Pe{X) = S(X — ^) as in the continuous case does the job. But this 
leads to zero probability for A^ not equal to an eigenvalue. We can correct for this by smearing out the delta function 
over a range of Xe so that Pe{Xe) is not zero between eigenvalues and a range of Xe corresponds to the same state. 
There are innumerable ways to parameterize & to achieve this. Here is one way. Define the & function 

= &,n n = n(Xe) (3.3) 

where Xe has no units and n(Xe) is the closest integer to X?. With this ^e(Xe) function we achieve agreement between 
Eq.(3.1) and Eq.2.1) using 

P t (Xe) = Pt(n(\t)). (3.4) 

where Pt(n(Xt)) = if there is no eigenvalue associated with the integer n(Xe). For discrete spectra with integer 
eigenvalues we can use the explicit forms 

oo 

Gi (X e (t)) = (n/(t) + 1/2 - X e (t)) P e (n e (t)) + ]T P e (j) 

j=n e (t) + l 
n e (t)-l 

U (X e (t)) = (X e (t) - n t (t) + 1/2) P t (n e (t)) + ]T ^0')> ( 3 - 5 ) 

j=-oo 

in the expressions for the current operator Jg. 

Using these definitions in the velocity expression Eq.(2.24) or Eq.(2.26) and keeping in mind that the physical 
values of ^ with discrete spectra are determined by Eq.(3.3) we have defined a deterministic Bohmian mechanics for 
operators with discrete and continuous spectra. The initial A configuration is taken from the quantum probability 
distribution Eq.(2.5) with Eq.(3.4) used for projectors for discrete operators. This means that the particular for a 
given ^£.„ is chosen at random from a distribution spread uniformly from n- 1/2 to n+ 1/2. We are free to order 
the eigenvalues along the Xe line any way we like, a freedom that also exists in the continuous case. This is the third 
way that the dynamics are not unique. The choice of ordering of the eigenvalues in A space for each beable operator 
profoundly effects the dynamics since, for example a system in eigenvalue state 2 can only get to eigenvalue state 4 
by first passing through eigenvalue state 3. For some systems the Hamiltonian determines a particular order, but 
perhaps there are Hamiltonians with transition elements between eigenvalues that are far apart in A space, for any 
ordering that one chooses. This is a possible disadvantage of this method which is not shared by the stochastic Bell 
scheme. 

In the next section we show how the formalism works for the simple example of Bohmian mechanics with only 
one beable and also present a visualization of Bohmian mechanics that allows us to dispense with the auxiliary Xe 
variables altogether. 
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IV. THE EXAMPLE OF ONE BEABLE AND A VISUALIZATION OF BOHMIAN MECHANICS 



Bohmian mechanics is integrable for the case in which there is only one operator promoted to beable status. For 
in that case 

< Q\dL(\(t),t)\0 >=d< 0|L(A(f),i)|0 >= 0. (4.1) 

so the solution is 

< 0|L(A(i), t)\0 >=< t\L(X(t))\t > = L Q . (4.2) 

The integration constant, L is uniformly distributed between zero and one. This equation has the following visual 
interpretation. Consider a line that goes from zero to one. The integration constant, L sits immovable on the 
line. Associate the portion of the line from zero to < t\P a \t > with beable value a, the portion of the line from 

< t\P a \t > to < t\P a \t > + < t\P b \t > with beable value 6, the portion of the line from < t\P a \t > + < t\P b \t > to 

< t\P a \t > + < t\P b \t > + < t\P c \t > with beable value c. These boundaries change with time. The value of the 
beable at any time is the value associated with the portion of the line that Lq sits on at that time. As an example 
consider a two state case in which an operator £, with eigenvalues ±1 and projectors P± is the beable (we have chosen 
a slightly different parameterization of the eigenvalues than we did in the previous section to take advantage of the 
symmetry in the two state case) . The Bohmian dynamics dictate that £ = —1 for < t\P-\t > > L and £ = +1 for 
io > <t\P-\t >. These results can be combined into the equation of motion 



or using < t\i\t >= 1 - 2 < t\P- \t > 



£(t)=sign(Lo-<t|P_|t>) (4.3) 



£(t) - sign(< t|£|t > -Co)- (4.4) 



where £o = 1 — 2Lq which is uniformly distributed from -1 to 1. Note that £o encodes observable information about 
the times tj that £ changes its state via < tj\^\tj >— £o so it is not "hidden", although it is uncontrollable. Also, the 
average value of £ (t) over all £o agrees with the quantum expectation value 



/: 



hl H(t)P(toWo=<t\i\t> (4.5) 



as is required for Bohmian mechanics to be consistent with quantum mechanics. 

The exact solution for one beable suggests a visual interpretation of Bohmian Mechanics for any number of beables 
that allows us to dispense with \g as we were able to do for the one beable case. For n beables consider an n 
dimensional space of area 1 with fixed boundaries. The space is divided into several n-dimensional bubbles. Each 
bubble corresponds to a particular E configuration and the volume of each bubble is the quantum probability of 
that configuration. Since the probabilities change with time, the bubbles are continuously contracting and expanding 
against each other. An immovable point is chosen at random in the n-dimensional space. The physical E configuration 
at time t is the E configuration corresponding to the bubble enclosing the immovable point at time t. 



V. SUMMARY 



In this paper we have generalized Bohmian mechanics so that it can incorporate beables associated with an arbitrary 
set of commuting continuous and discrete operators. The equations are deterministic and time reversible and agree 
with Bohm's original formulation for the case of continuous position operators. The simple case of only one beable is 
presented and the solution suggests an intriguing visualization of Bohmian mechanics. 
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